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Abstract. We study a one-dimensional extended Hubbard model with longer-range Coulomb interactions 
at quarter-filling in the strong coupling limit. We find two different charge-ordered (CO) ground states as 
the strength of the longer range interactions is varied. At lower energies, these CO states drive two different 
spin-ordered ground states. A variety of response functions computed here bear a remarkable resemblance 
to recent experimental observations for organic TMTSF systems, and so we propose that these systems are 
proximate to a QCP associated with T — charge order. For a ladder system relevant to Sri4Cu24 04i, 
we find in-chain CO, rung-dimer, and orbital antiferromagnetic ordered phases with varying interchain 
couplings and superconductivity with hole-doping. RPA studies of many chains (ladders) coupled reveal 
a phase diagram with the ordered phase extended to finite temperatures and a phase boundary ending 
at a quantum critical point (QCP). Critical quantum fiuctuations at the QCP are found to enhance the 
transverse dispersion, leading to a dimensional crossover and a T = decofinement transition. 

PACS. 71.30.-|-h Metal-insulator transitions and other electronic transitions - 71.28.-|-d Narrow-band 
systems - 72.10.-d Theory of electronic transport 



' 1 Introduction 

Electron crystallization, or the charge ordering of elec- 
trons due to interactions, is an issue of enduring inter- 
, est in condensed matter physics [l^- A host of excel- 
' lent studies clearly show the relevance of charge ordering 
I (CO) in diverse systems like organics [3], transition metal 
. oxides mSj, and coupled chain-ladder systems [S]. Gener- 

■ ically, CO seems to compete with superconductivity [3l[5], 
, or with metallicity OH]. Hence, in addition to its intrin- 

■ sic academic interest, the study of the conditions favoring 
I CO, along with its competition with metallic (magnetic) 

and/or superconducting states constitutes a problem of 
wide-ranging interest for a host of real systems. 

The simplest case of a CO state is the one where band 
fermions hop on a lattice with a staggered (static) crystal 
potential [6j. In this case, the CO gap is just the band- 
insulating gap. Needless to say, this picture is too sim- 
ple; it does not contain the ingredients necessary to study 
the coupled charge/spin aspects of the systems mentioned 
above: these point to the basic importance of the Hub- 
bard U in practice. In contrast to the band-CO case, one 
expects the correlation driven CO state to have drastic 
consequences for the magnetic order emerging at lower 
energy scales [2!i!3l'4l'5] , as well as for the competition be- 
tween charge and spin orders manifested as one between 
Mott insulating (magnetic) and metallic (superconduct- 
ing) states. 



The availability of several powerful analytical as well 
as numerical methods in one dimension facilitates the de- 
tailed study of the effects of strong correlations in such 
systems. The study of the Hubbard model in one dimen- 
sion, with and without extended Coulomb interactions, 
has a long history [TIH] • With the focus mostly on the case 
of 1/2-filling [9], systems away from 1/2-filling have not 
received sufficient attention. Furthermore, studies at dif- 
ferent fillings have typically concentrated on the weak cou- 
pling limit [igifTTlfT^ . while real systems of interest 
[S] are generically in the strong coupling limit of appropri- 
ate Hubbard-type models. Additionally, one dimensional 
correlated models with couplings to the lattice have also 
been studied [T51fTH[T5] . revealing phase diagrams with 
various different kinds of charge orders (CO) and spin or- 
ders (SO). 

The issue of whether various types of CO can, at fill- 
ings away from 1/2, be driven purely by longer range 
Coulomb interactions remains, however, largely unadressed. 
The importance of extended Coulomb interactions was 
appreciated in an early work by Hubbard [16] . Here, it 
was shown that the next nearest neighbor Coulomb (nnn) 
interaction can be as large as 40 percent of the nearest 
neighbor (nn) interaction, with the intra-atomic Hubbard 
U being the largest, and the one-particle hopping between 
neighbors the smallest energy scales in an effective one- 
band extended Hubbard model. In practice, coupling to 
high-energy Einstein (optical) phonons has the result of 
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reducing the nn coupling |14j ; in a model with Coulomb in- 
teractions extended to next nearest neighbours, this gives 
rise to the possibility of tuning the ratio of nnn to nn in- 
teractions through 1/2, a point at which interesting addi- 
tional physics is expected. Thus, the search for competing 
charge ordered states driven purely by long-range elec- 
tronic interactions, and the co-existence of charge and spin 
order, in a one dimensional 1/4-filled electronic model are 
among the primary goals of the present work. Further, by 
working with a model in which the Coulomb correlations 
are much larger than the single particle hopping strength, 
we are able to obtain analytic results in a regime where 
little progress has been made. From our earlier discussion, 
our results are also clearly relevant to several material sys- 
tems. 

In the following section, we begin by presenting a deriva- 
tion of a ID transverse field Ising model (TFIM) effective 
pscudospin Hamiltonian in the extreme strong-coupling 
limit for the charge degrees of freedom of the 1/4-filled 
ID spinful electron Hubbard model with extended corre- 
lations. We then discuss the consequences various features 
of the phase diagram of the (exactly solvable) TFIM effec- 
tive model have for the charge and spin degrees of freedom 
of the original electronic model. Various thermodynamic 
quantities of the electronic model, e.g., the optical con- 
ductivity, the dielectric function and the electronic contri- 
bution to the Raman scattering, are computed from the 
effective pseudospin theory. We follow this, in section III, 
with an analysis of the effective theories obtained for the 
case of a 2-leg ladder of two coupled TFIM chains with 
the inter-chain couplings in the strong coupling limit, as 
well as for a ladder doped with holes. In section IV, we 
then analyse the case of many such TFIM (chain and 2- 
leg ladder) systems coupled to another by a random phase 
approximation (RPA) method [T UTfUTglfTO] . The goal is to 
conduct a systematic study of the physics of dimensional 
crossover and deconfinement in such anisotropic strongly 
coupled lower dimensional systems, when coupled to one 
another. In section V, we present a comparison of our 
findings with some recent numerical works. Finally, we 
conclude in section VI. 



2 Single chain in strong-coupling regime 

In this work, we study this issue within an extended quarter- 
filled Hubbard model on a linear chain, described by, 

i/cff = -tY,iCla+i,^ + h.c) + iU- 2zP)J2n,^n,i 

i.a i 

+V^nini+i + (1) 

i i 

The couplings U , V and P represent the on-site Hubbard, 
nn and nnn replusive interactions, t represents the hop- 
ping strength and z the coordination number for the lat- 
tice (which in the case of the linear chain is 2). This some- 
what complicated looking microscopic model gives us the 
means to model extended Coulomb interactions in a 1/4- 
filled system in the limit of strong coupling (as discussed 



below). Further, it is by now a well accepted paradigm 
that in theories of electrons in one-dimension interacting 
with one another through well-behaved (i.e., non-singular) 
potentials, the low-energy effective (fixed point) theory is 
one where the charge and spin degrees of freedom have sep- 
arate dynamics [7] . This is the concept of the Tomonaga- 
Luttinger liquid with spin-charge separation, and charac- 
terised by the spin fluctuations being those of an ideal 
5"= 1/2 XXX AF chain 

77^ Sj .Sj+i, (2) 

while the charge fluctuations are described by the Hamil- 
tonian 

He = -t^(4c,+i + h.c) + {V - J 1 4) mn.+i 

i i 

(3) 

i 

that describes a model with frustrating interactions. The 
coupling J is obtained from straightforward perturbation 
theory as J 4t^/(C/ - V - AP). In ID, the projected 
fermions are spinless fermions with a hard-core constraint. 

2.1 Effective TFIM model 

The model for spinless electrons given above has been con- 
sidered as a model for studying the effects of frustration 
on charge ordering |201l21j . For narrow-band systems, we 
consider the limit t << {V — J /A), P. Before embarking 
on a detailed derivation of the effective Hamiltonian in 
this regime, we present a heuristic one [20]. For this, we 
employ an extension of the trick used for the Id next- 
nearest neighbor Ising chain: for t = 0, we notice that 
with {V - J/4) > 2P, the ground state is the usual CDW 
(Wigner) crystal (and doubly degenerate, (..101010...) and 
(...010101...), for 1/4-filling). With 2P > (V- J/4), how- 
ever, the dimerized state (Peierls) is the ground state. 
The ground state has a degeneracy of four, with one of 
the states written schematically as (11001100....) and the 
other three achieved by shifting the state by one lattice 
site. Taking the state shown and splitting this in a slightly 
different way, we have [...(01)(10)(01)(10)...]. Associating 
a pseudospin r = 1/2 operator, with = +1 for (10) 
and -1 for (01), the state is antiferromagnetic and doubly 
degenerate in terms of the . This state has a partner, 
obtained by shifting the state given above by 2 sites (or 
the pseudospin r variables by one site). For small t, this 
is an attractive trick because (in spin language) the trans- 
verse term does flip the rf ^ but cannot break a pair. This 
leads to the effective Hamiltonian p01l21| . 

iJeff = - E [2^^r + (V^ - ^/4 - 2P)rf rf+i] (4) 
I 

Let us now see how this simple Hamiltonian is obtained 
as the effective Hamiltonian for our original theory for the 
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charge degrees of freedom in the hmit of the couplings 
U >> V, P >> t. We begin with the Hamiltonian for a 
spin chain system 

n 

+J2S^nSU^ + hS^^] (5) 

where S^, and and spin- 1/2 operators. The cou- 
plings Jx, Ji, J2 > are the nearest neighbour (nn) XY, 
the nearest neighbour Ising and the next nearest neigh- 
bour (nnn) Ising couphngs respectively and h is the exter- 
nal magnetic field. For h = 0, this Hamiltonian can be de- 
rived from the on-site Hubbard interaction U —^ 00 strong- 
coupling limit of an extended Hubbard model at 1/ 4-filling 
and with nn ( Ji = {V — J/4)) and nnn (J2 = 2P) density- 
density interaction couplings and the nn electron hopping 
(Jx = t) via a Jordan- Wigner transformation (from a 
model of spinless-fermions to spins) 

H ^^[-^{clci+i+h.c) + Jinin,+i + J2nin,+2] ■ (6) 

n 

We study the problem in the limit of strong-coupling where 
Ji, J2 » Jx (but where (Ji — 2J2) ~ 2Jx)- 

Let us begin by studying the case of J^. = (we 
will be studying eq.© for the case of ft. = in all that 
follows). It is easy to see that for the case of Ji > 2J2, 
the ground state of the system is given by a Neel-ordered 
antiferromagnetic (AF) state with two degenerate ground- 
states given by 

\AFGS1) = 1... + - + - + - + -...) 
\AFGS2) = I ...- + - + - + - + .. .) (7) 

where we signify = 1/2,— 1/2 by + and — respectively 
and we have explicitly shown the spin configuration in the 
site nos. —3 < rt < 4 in the ground states. In the original 
electronic Hamiltonian eq.([6]), this AF order corresponds 
to a Wigner charge-ordering (CO) in the ground-state. 
Similarly, for the case of Ji < 2 J2, the ground state of the 
system is given by a dimer-ordered (2, 2) state [20] with 
four degenerate ground-states given by 

I22G51) = 1...- + + -- + + -...) 
\22GS2) = 1...-- + + -- + + ...) 
I22G53) = 1... + -- + + -- + ...) 
I22G54) = 1... + + -- + + -- ...) 

(8) 

where we signify = 1/2, —1/2 by -I- and — respectively 
and we have explicitly shown the spin configuration in the 
site nos. —2 < n < 5 in the ground states. In the original 
electronic Hamiltonian eq.®, this (2, 2) order corresponds 
to a Peierls CO in the ground-state. Further, as the Hamil- 
tonian eq.([5]) with Jx = = h is the one-dimensional axial 
next nearest neighbour Ising (ANNNI) model, we know 
that the "frustration" point at Ji — 2J2 [12] is gapless, 
highly degenerate and has a ground state entropy per site 



of ln(l + \/5)/2 ~ 0.4812 while the ordered ground states 
on either side of the frustration point have zero ground 
state entropy |23j. 

We now work out the effect of the XY terms in the 
Hamiltonian (O on these ground states. Let us start with 
noting the effect of a XY term on a single nn spin-pair 
on the 4 degenerate (2, 2) ground-states; for purposes of 
brevity, we will denote the entire Jx{S!^S^_^_i + Sn^S^j^-^) 
term simply as J"'"+^. Thus, 



J°'i|22G51) 


= 4^'\.- 


. + -...)^||. 


> 


J°'^|22G52) 


= J'A-. 


. + + ...) =0 




J°'i|22G53) 




.- + ...)-||... + -.. 


) 


J°'i|22G51) 


= J°A-- 


.--...)=0 




Jx^'^\22GSl) 




.. + + ...) =0 




J^^'°\22GS2) 






• •> 


4~i'"|22G53) 




..--...)=0 




Jx^'°\22GSA) 






..)(9) 



In a similar manner, we study the action of the opera- 
tor J"'""*"^ on the two degenerate ground states of the AF 
ordered configuration as 



Jl^^\AFGS\) 


^J'x'\... + -. 




+ ...) 


Jl^^\AFGS2) 






-...) 


Jx^'^\AFGS\) 


= Jx'^'\...- + 


' 2 ' 


f-...) 


J-i'°|AFG52) 




...)= ^|...- 

' 2 ' 


- + ..01O) 



Defining bond-pseudospins t/ = [Sf — S^_i)/2, r+ = 
S^S~_i and t~ — S~ S^_j^ (which can be rewritten in 
terms of fermionic operators in the original electronic Hamil- 
tonian eq.® as rf = (n^ — ni_i)/2, = c\ci-i and 

= Cicl_-^ respectively), we can write the four degen- 
erate ground states of the (2, 2) ordered configuration in 
terms of these bond pseudospins as 

|22GS'l) = |...0- 0-1-0...) 

|22GS'2) = |... + 0- 0+ ...) 

|22GS'3) = |...0 + 0-0...) 

|22GS'4) = |...-0 + 0- ...), (11) 

where we have denoted — 1/2 as + and — —1/2 
as — and have explicitly shown the pseudospin configu- 
rations on the bond nos. < n < 4. We can clearly see 
from eg. pTj) that these four ground-states break up into 
two pairs of doubly degenerate (AF) orderings of the pseu- 
dospins defined on the odd bonds (|22G51) and |22GS'3)) 
on the even bonds (|22G52) and |22GS'4)) respectively. 
It is also simple to see from eq.® that the action of the 
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operator J"~^'" (for the nearest neighbour pair of sites 
given by (n — l,n)) on these four ground states is to flip 
a pseudospin defined on the bond n (lying in between the 
pair of sites (n — 1, n)) or to have no effect at all. 

We can now similarly see that the two-degenerate ground 
states of the AF ordered configuration can be written in 
terms of the bond-pseudospins defined above as 



\AFGSl) = 
\AFGS2) = 



(12) 



where we have explicitly shown the pseudospin configura- 
tions on the bond nos. < n < 4. From eg. (fT2| . we see 
that the two degenerate ground-states have antiferromag- 
netic ordering of pseudospins on nn bonds; this can equally 
well be understood in terms of the ferromagnetic ordering 
of pseudospins on the odd bonds and on the even bonds 
separately. Further, from eq. lTUl) . we can see that the ac- 
tion of the operator J""'^'" (for the nearest neighbour pair 
of sites given by (n — 1, n)) on these two ground states is 
again to flip a pseudospin deflned on the odd (even) bond 
n (lying in between the pair of sites {n — l,n)) against a 
background of ferromagnetically ordered configuration of 
pseudospins defined on the odd (even) bonds. 

Thus, we can model these pseudospin ordered ground 
states (|11I) . (|T2|) as well as all possible pseudospin-flip ex- 
citations above them (as given by action of operators of 
the type J"~^'" ®,(Iini)) with the effective Hamiltonian 

H^-J2 + (^1 - 2 J2)«+2] 

n^odd 



n+2J 



. [2J.< + (Ji - 2J2)<r, 

^[2J.< + (Jl-2J2)«+2] , 



is qualitatively similar to that of the usual t — J model. 
The "Mott" insulating state for F - J/4 > 2t has Wigner 
CO in the ground state, and the M-I transition is of the 
Kosterlitz-Thouless type [29] . 



T 



Low T gapped 
CO state 
with LRO 



Quantum 
Critical 



Low T gapped 
short ranged 
charge VB liquid 

9 ' 



(13) 



Fig. 1. A schematic phase diagram of the temperature T vs. 
coupling 5= |(V- J/4-2P-2t)/(V- J/4-2P)| = Zi,/(V- 
J/4 — 2P) of the effective ID Transverse Field Ising Model 
theory for the charge sector of our original electronic model. 
The black circle with g = Q represents the T = Q quantum 
critical point (QCP) separating an phase with true LRO given 
by Wigner CDW or Peierls dimer order (depending on whether 
V-J/4: > 2P or V^- J/4 < 2P respectively, thick dark line) and 
a quantum disordered phase corresponding to a charge valence 
bond (VB) liquid. The finite-T physics of the quantum critical 
region lying just above the QCP is discussed in subsection 
IIB. The dashed regions represent finite-T crossovers to low-T 
gapped charge ordered and charge VB liquid phases with no 
long-range order (LRO). 



where n is the bond index. 

This is just the Ising model in a transverse held, which 
is exactly solvable [24l[25] and has been studied extensively 
in ID [IZllin]- If {V - J/4 - 2P) > 0, the ground state 
is ferromagnetically ordered in r^, i.e, it corresponds to a 
Wigner CDW. For {V - J/4 - 2P) < 0, the Peierls dimer 
order results in the ground state. At {V — J/A — 2P) < 
2t, the quantum disordered phase has short-ranged pseu- 
dospin correlations, and is a charge "valence-bond" liquid. 
The quantum critical point at {V — J/A — 2P) = 2t sep- 
arating these phases is a deconflned phase with gapless 
pseudospin (r) excitations, and power-law fall-off in the 
pseudospin-pseudospin correlation functions. Correspond- 
ingly, the density-density correlation function has a power- 
law singular behavior at low energy, with an exponent 
a = 1/4 characteristic of the 2D Ising model at critical- 
ity. The gap in the pseudospin spectrum on either side of 
the critical point is given hy Ar = 2\V - J/4 -2P -2t\. 
Further, the quantum-critical behaviour extends to tem- 
peratures as high as T ~ Ar/2 [28] and undergoes flnite- 
temperature crossovers to the two gapped phases at T ~ 
\Ar\. For P = 0, the metalhc phase for V - J/4 < 2Hs a 
Luttinger liquid, and in this limit, the low-energy physics 



The full Hamiltonian in our case for the strong-coupling 
limit is now given by 

iJeff =-^[2irf+(y-J/4-2P)Tfrf+i]+J^S,.S,+i (14) 



To study the magnetic phases, we adapt the Ogata- 
Shiba [30] technique for our case. This is possible if J(~ 
/[U — V — 4P)) << t, V, in which case, the pseudospin 
part is first solved exactly (this is possible because of 
the known exact solution of the ID transverse field Ising 
model), and the exchange part is then treated as a pertur- 
bation. Writing the total wavefunction as a product of a 
spin and pseudospin wavefunction (where the spin wave- 
function is defined, for a system of size N, in a Hilbert 
space of dimension 2^), i.e, \tp >= |r > (>^\S >, and fol- 
lowing standard degenerate perturbation theory, the spin 
degeneracy is lifted by the correction (of order l/N): 



< iJ^ff >'= -2t < > 



I 



Ju+i(Si.Si+i-l/4) (15) 



where the average < .. >' denotes that the average is taken 
over the exact ground state |r > of the pseudospin part 
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hiv- 



N 



above, i.e. < A >'—< t\A\t > and Ji ;+i 
2P) < rfrf+i >. 

An interesting fact now emerges: Wigner CO (FM or- 
der of r) results in an HAFM S = 1/2 spin model with 
the Hamiltonian Hg — S^.S^^-i. This gives rise to 

a gapless AF ground state for the spin degrees of free- 
dom. The charge (holon) excitations have a gap At = 
2{V — J/4 — 2P — 2t); this corresponds to a linear confin- 
ing potential for holons. On the other hand, Peierls dimer- 
ization in the charge sector (AF Neel order of r) gives 
rise to dimerization in the spin sector, with the Hamil- 
tonian = JJ2t[^ + (-l)'(5]S,.S,+i. The bare value of 
the dimerisation is given found from the perturbation the- 
ory as proportional to the Ising pseudospin (charge dimer) 
gap, i.e., S ^ Ar/J- However, as we shall discuss now, a 
continuum field theoretic description of this Hamiltonian 
shows that the dimerisation parameter S is renormalised 
by quantum fluctuations, and scales as (S ~ (Z\t-/ J)^^"^. 

Translated into fermion variables, this yields a sine- 
Gordon problem with = 27r, and describes an instabil- 
ity to a singlet pinned ground state commensurate with 
the Peierls CO setting in at higher energies. The elemen- 
tary excitations are solitons carrying = ±1. Scaling 
theory predicts a dimer gap, Ar ~ S^^'^J. Exactly at 
P"^ = 2n, the SG model has just two S"^ = breather 
excitations with opposite parity |31| . the lowest, even par- 
ity breather being degenerate with the = ±1 soliton 
doublet, forming a S" = 1 triplet, while the second odd- 
parity breather is a singlet with a gap, \/3Ar. It is impor- 
tant to notice that both charge and spin order arise from 
long range Coulomb interactions, and do not involve an 
electron phonon coupling mechanism. 



2.2 Consequences of CO on thermodynamic quantities 
in the quantum critical regime 

Let us consider the implications of having the CO state in 
the high-T regime, where one could imagine the system to 
be effectively one-dimensional. In particular, we want to 
look at the to, T dependence of the various response func- 
tions at high-T. Using the exact solution of the pseudospin 
model in ID, the high T (in the "quantum critical" regime) 
behavior can be explicitly derived [27 . In fact, near Ising 
criticality, the response function, xi^) — r^^/^ where r — 
{x^ -^T^Y^"^ (with the velocity v set to unity). This rela- 
tion is still valid away from criticality in the "short range" 
region, r « Ar, where Ar = 2\V - J/4 - 2P - 2t\ is 
the pseudospini gap (which corresponds, in our case, to 
the charge gap) of the ID-TFIM. As we will see, in what 
follows, knowing the microscopic form of the Ising pseu- 
dospin gap At- is enough information for us to be able to 
compute a host of thermodynamic response functions of 
the sytem and be able to relate them to the couplings in 
our original microscopic theory of interacting electrons in 
one dimension. Using this asymptotic form, we have 



Xcrit(0,w) 



{2nT) 



(16) 



where S — and A — 1/16 is the conformal dimension. 
B{x,y) is the beta function. 

In the quantum critical region, an illuminating form is 



x(o,o) 



(17) 



where the quantities Tr, = (2 tan(7r/16)fcBr/^)e~'^"/''^'^, 
wi = 0.795(fcBT/;i) and ^ = 1.28(cft/fcBr)e^-/'=^''^, are 
determined solely by T and the fundamental natural con- 
stants, as expected in the QC regime. We stress once again 
that here, Ar = 2(t/- J/4 - 2P- 2t) is the energy gap to 
charge excitations in the Wigner/Peierls CO states (Ising 
pseudospin gap) described above. This represents the col- 
lective charge susceptibility, and the optical conductivity 
follows directly from a-{uj) = —iuxiO^w), giving. 



x(o,o) 



(18) 



The frequency dependent dielectric function, e{uj), is 
obtained from e(w) = 1 + {4:TTia{uj)/uj), and the electronic 
contribution to the Raman scattering is estimated there- 
from to be given by I-r,{u!) ~ lm(l/e(0, w)), for light po- 
larized along the chain axis. In terms of the charge sus- 
ceptibility, this is simply. 



/r(ci;) = Im- 



1 



(19) 



where F{uj,T) = jr. , , ^ ^^ , . — , „ .o . 

x" {k, tjj) has its maximum value at ojm = wi— i(w^//r,), 
implying that the collective mode broadens and shifts to 
higher energy linearly in T with increasing T at high tem- 
peratures. Further, the T-dependent damping rate of the 
collective mode correlates well with the relaxational peak 
seen in transport, underlying their common origin. In fact, 
the dc resistivity is linear in T at high T, with "insulat- 
ing" features showing up at lower T . In our picture, these 
are collective (longitudinal) bosonic charge-density modes 
in the high-T quantum critical region above an incipient 
CO transition (expected to occur at low T). In fig.(l), we 
show the electronic Raman lineshape as a function of w/T. 
The sharp low energy peak corresponds to the collective 
charge density fluctuation mode of the CO ground state. 

The corresponding frequency-dependent dielectric con- 
stant also shows an explicit oj /T scaling in the QC regime, 
or generally, at high-T, it shows strong T-dependence. 
From fig. (2), we see that it becomes w-independent at 
high T, but appreciably increases as T is lowered, with 
a maximum at w ~ T. 

The fact that organic charge transfer salts [22i exhibit 
features very similar to those found above has interest- 
ing implications. In light of our results, these anomalous 
features can now be identified with proximity to an un- 
derlying quantum critical point associated with charge 
(Wigner/Peierls) ordering. We recall that very recent work [33l 
IM] shows that the dimerized insulating state in TMTSF 
systems has charge order at low T. Interestingly, e'(0,w) 
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CO/T 



Fig. 2. A three-dimensional plot of the Raman intensity 
In versus a scaled energy ui/T and temperature T for pa- 
rameter values of the original model which constitute a gap 
= 0.05fcB. 



a((B)/Q) 




m/T 



Fig. 3. A three-dimensional plot of a{uj)/uj versus a scaled 
energy uj/T and temperature T for At = COSfce- The limiting 
behaviors are: (1) o-{lo) ^ jItt74 high T, and a{u)) ~ ^-^^ for 
low T. (2) e'{Lu,T) = 1 -I- ^ for high T, and e'{co,T)^~ ^ 
for low T. 



indeed shows appreciable increase as T is lowered, further 
supporting an interpretation based on proximity to an un- 
derlying CO ground state [32] ■ Hence, we conclude that 
observation of these features in TMTSF systems consti- 
tutes strong evidence that this system is close to a putative 
QCP associated with charge order. Observation of dimer- 
ized/Neel ordered AFM states co-existing CO states at 
low T are also naturally understood in light of the analy- 
sis above (341. 



3 Two-chain TFIM Ladder systems at strong 
inter-chain coupling 

We now consider the strong coupling version of a coupled 
two-chain ladder system, with each chain being described 
by H as in eq.(l). In the strong coupling limit, where each 
chain is described by a TFIM for charge degrees of free- 
dom, the coupled chain model is constructed as follows. 



For C/ ^ cx), and y, F > i (but {V ~ J/4 - 2P) compa- 
rable to t) , the charge degrees of freedom of the fermionic 
problem for each chain are described by an effective pseu- 
dospin model on n-n bonds, via the effective Hamiltonian, 



j 



[2tT^ + {V- J/4 - 



2P) r]r]^^ 



(20) 



Rotating the pseudospin axis such that t^, — 

— and coupling two such chains via an interaction cou- 
pling Ui_ and a two-electron interchain transfer we 
have the effective Hamiltonian for the charge sector of the 
two chain system as 



H 



{V~J/A-2P)tI,T 



,)] (21) 



where a, 6 = 1,2 is the chain index. Denote the in-chain 
pseudospin coupling as J = (y — J/4 — 2P) and the inter- 
chain pseudospin coupling as J± = U±. Here, we study 
the strong coupling version of this problem in several lim- 
its by deriving the respective low-energy effective Hamil- 
tonians (LEEH). The weak-coupling problem is studied 
elsewhere l35l. 



3.1 The case of |J±| >> |J|, i± 

For the case of | Jj_| >> |J|, t± , the 2 chain system can 
be better thought of as strongly-coupled rungs which are 
weakly coupled to their neighboring rungs. Thus, we treat 
J as a weak perturbation on the zeroth-order system of 
rungs defined by the large coupling J±, giving iJe// = 
Hq + Hi where 



Ho 
Hi 



J 



3,a 



5Z '''3,a'''j,b 



h.c) (22) 



j.a j.a.b^a 

where the effective magnetic field is given by /i = 2i > 0. 
3.1.1 LEEH for J_l < 

For Jx < and h « J±, wc find that the triplet state 
1+) = 7f(l Ti) + I it)) and the singlet state |-) = ^(| Ti 
) ^ I it)) are degenerate on any rung and are separated 
from all other states by a large gap of order J±. Thus, 
these two states define the subspace which will deter- 
mine the low-energy physics of the system. Identifying a 
pseudospin- 1/2 operator with the low-energy subspace 
on each rung, we treat the Hamiltonian Hi as a pertur- 
bation (to second order in J/J±) to obtain the LEEH as 



H 



2J, 



E 



2J, V? 



2h^ 



Ah- 



t± 



3 



0-1. ^J2 -4/l2^^^^' 



72 t2 



2h^ 



4/i2 



)] (23) 
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We find that tj_ acts as the strength of a Zeeman-splitting 
like term in the LEEH. Bosonising this, we obtain a sine- 
Gordon Hamiltonian with a cosine potential in the dual 
(6) field and a magnetic-field term 

H = ^[{d.^f + {d.ef] - ^ cosM - ^a,</> , (24) 
2 Zna 2 

where the cosine term arises directly from the BCS-like in 
the effective Hamiltonian, with the bare value of the mass 
m ~ J^h? / J^. The velocity, u, and sine-Gordon correla- 
tion exponent, (3i, are both functions of the energy scale 
J'^/J_L. We note that bosonisation of the general XYZ 
Hamiltonian results in the appearance of an additional 
Akp Umklapp term [7], cos/32'/', which is irrelevant for a 
finite t± and is hence ignored in what follows. When t±_ 
is below a certain critical value, incommensurate Wigner 
charge order (ordering of the ^^) occurs |36]. Above this 
critical value, a spin-flop transition orders the system in 
the X direction (i.e., ordering of the ^^) via a Kosterlitz- 
Thouless transition. For 01 < Stt, the cosine in the dual 
field, 6, is a relevant perturbation and orders the dual field. 
The magnetic-field term cx t± leads to a ground state with 
charges which are coherently delocalised on the diagonals 
of each pair of nearest-neighbor rungs; this is an orbital 
antiferromagnet-type ordering with circulating currents in 
plaquettes [29ll7] . 

3.1.2 LEEH for Jj_ > 

For J± > 0, and h > 0, the triplet state |+) = | tt) is 
the low energy state on any rung. For /i = 0, we find that 
the triplet states |-|-) (defined above) and |— ) = | ii) 
are degenerate. Thus, we can again identify these two 
states as the subspace which determines the low-energy 
physics of the system. For h « J, we again identify a 
pseudospin-1/2 operator with the low-energy subspace 
on each rung, and treat the Hamiltonian Hi as a pertur- 
bation (to second order in J/J±) to obtain the LEEH as 

^ = -ljE^Ki+i-2tE^|- (25) 

i 3 

This is just the ID TFIM (with ferromagnetic Ising cou- 
phng). In the ordered phase, the ground state has in-chain 
Wigner CO and dimers on every alternate rung. The dis- 
ordered phase is a gapped, short-ranged charge-dimer liq- 
uid. At t = J'^/AJ±, the quantum critical point describes 
a gapless charge-dimer liquid with u/T^vk/T QC scal- 
ing, exactly as was described before. Transposing the re- 
sults obtained before, we conclude that the dc resistiv- 
ity, optical conductivity, electronic Raman and dielectric 
responses will be exactly described by the same scaling 
functions (eqs. p^ - fTOl) ) with the gap, Ar, now being the 
CO gap of the ladder problem {H in eq. p3|) '). Very in- 
terestingly, exactly such behavior is observed in undoped 
ladder system SriiCu2iOii [37] and attributed to a lon- 
gitudinal, collective charge fluctuation mode, exactly as 
described here. 



3.2 \t^\ » |J|, 

For the case of | Jj_| >> | J|, ij^, we can again treat the 2- 
chain system as that composed of strongly coupled rungs 
which are weakly coupled to their neighbours. Thus, we 
treat J as a weak perturbation on the zeroth-order system 
of rungs defined by the large coupling Jj^ , giving H^f / — 
Hq + Hi where 

Hi = -jT.^lar-+i,a-J± E ^larlb (26) 

j,a j.a.b^a 

where the effective magnetic field is given by /i = 2i > 0. 

For tj_ > and h = t±/2, we find that the triplet zero 
state |-|-) = ^(1 ti> + l it)) and the triplet up state |-) = 
I tt) are degenerate on any rung and are separated from 
all other states by a large gap of order t±. Thus, these two 
states define the subspace which will determine the low- 
energy physics of the system. Identifying a pseudospin-1/2 
operator with the low-energy subspace on each rung, we 
treat the Hamiltonian Hi as a perturbation to obtain the 
LEEH as 

H - Et-'^^^l + f (^/^7+i + h-c)] (27) 

3 

Note that, unlike the LEEHs derived earlier, this LEEH is 
at first order in J/t± and J±/t±. Further, we have checked 
that the terms obtained at next order in the perturbative 
expansion are considerably smaller and do not introduce 
anything qualitatively new; it is, therefore, sufficient to 
stop at this order. 

The expression (l?f)) is the Hamiltonian for the isotropic 
XY model in a transverse magnetic field. This theory is, 
again, exactly solvable and has a T = QCP in its phase 
diagram at = |J|/2. Further, there exists an equiv- 
alence between the classical 2D Ising model at finite T 
and the T = isotropic quantum XY model in a trans- 
verse field [35] ■ the high (low) T regions of the former with 
T > Tc {T < Tc) map onto the high (low) J regions of the 
latter with J± > {J± < J^). The critical exponents 
associated with the finite-T thermal phase transition in 
the classical 2D Ising model are also identical to those as- 
sociated with the T = quantum phase transition in the 
XY model in a transverse field. Thus, we can expect scal- 
ing forms for the various response functions of the system 
in the quantum critical region of the phase diagram (lying 
just above the QCP in the T direction), as discussed ear- 
lier. In the ordered phase, the ground state of the system 
has either rung-dimer or rung-diagonal dimer order, while 
in the disordered state is characterised again by a gapped 
charge-dimer liquid. 

3.3 LEEH for hole-doped ladder 

Upon doping the ladder with holes, while a single hole ex- 
periences a linear confining potential in the Wigner (Ising- 
like) or Peierls (dimerized) CO background, a pair of holes 



8 



Siddhartha Lai, Mukul S. Laad: Ordering from frustration in a strongly correlated one-dimensional system 



on the same rung is free to propagate. One can then de- 
scribe the hole-pair as a hard-core boson, representing its 
creation and annihilation operators using the spin-1/2 op- 
erators (T^; the local charge density is then described by 
(T^. Following [5^, we find the LEEH describing the dy- 
namics of such hole-pairs to be the XXZ model in an ex- 
ternal magnetic field 



h.c) - MhCTlcrl+i - AiCr|] (28) 



where ih ~ J'^/Ji^ is the pair- hopping matrix element, 
Uh is the Coulomb interaction between pairs on nearest- 
neighbour rungs and is the chemical potential of the 
holes. The phase diagram of this model is known [29]: for 
/Li = and Uh > ih, the ground state is an insulating 
CDW of hole pairs. Beyond a critical /ic = /(wh,ih), the 
system has a ground state described by Bose condensa- 
tion of hole pairs. In fact, from the bosonisation analy- 
sis of the equivalent S — 1/2 XXZ model in an exter- 
nal Zeeman field, we know that < crfCTf, ^ >~ and 



" where a = 1/2 



(2uh/th). 



Clearly, for a < 1, the ground state has dominant su- 
perconducting correlations. This is true for both the cases 
described above: in the first case, we have a Bose conden- 
sate of intrachain pairs of holes, while in the second hole 
pairs on individual rungs Bose condense, describing two 
possible superconducting types in the ladder system. This 
finding matches our conclusions obtained from a weak cou- 
pling analysis |35| , and thus constitutes a generic feature 
of undoped/doped strongly correlated ladder systems. 



is simple (see Fig.([T|)), with an ordered phase for 2t < 

V — J/4:—2P (g < 0), a quantum disordered phase for 2t > 

V — J/4 — 2P {g > 0) and a quantum critical point at 2t = 
V—J/A—2P {g — 0). A finite transverse couphng, denoted 
by 0±, causes the ordered phase to be extended to finite 
temperatures (with a critical temperature Tc for the case 
of 2t = V — J/A — 2P) with a first order phase boundary 
ending at a new quantum critical point (QCP) gc — '^c/V 
at T = 19]. As for the simple TFIM 27], there exists 
a "quantum critical" region just above the QCP and to 
the right of the ordered phase, with a crossover to the 
disordered phase at finite T. This is shown schematically 
in the T — g phase diagram given below in Fig.(j4|). The 
transition belongs to the 3D Ising universality class while 
the QCP to the 4D Ising universality class [39] . 




Ciiiplcss /' //' 
(Jiiniiliiiii /',' 
Critical ,' 

Ilegipn,-/' Disordered 
/' / .' Phase 



gc (QCP) 



4 Dimensional Crossover in Coupled Chain 
and Ladder Systems 

Having studied the case of the single chain at strong- 
coupling as well as those of strongly coupled 2-chain lad- 
der systems of such chains, we now turn our attention 
to the invesigation of the case of when many such chains 
and ladders are coupled to one another through transverse 
couplings of the kind studied in the previous section. We 
have recently studied [35] in detail the consequences of 
transverse couplings between many TFIM systems using 
the random phase appromximation (RPA) [T7[fT8lfT9] . In 
what follows, therefore, we will rely largely on the deriva- 
tions and results give there [32] . With the transverse field 
Ising model (TFIM) having appeared as an effective the- 
ory for both the case of a single electronic chain at strong 
coupling as well as proving to be the LEEH for a strongly 
coupled 2-leg TFIM ladder system, we will focus first on 
dealing with the case of TFIM chains and ladders. We 
then pass to a study of some of the other ladder LEEH 
theories for another work. 

Let us begin by setting out what we can expect to 
happen generically upon coupling many TFIM chain or 
ladder systems. In the absence of any transverse couplings, 
as discussed earlier, by defining the coupling g — \{V — 
J/4-2P-2t)/(F-J/4-2P)|,the phase diagram at T = 



Fig. 4. The T — g phase diagram for the case of many TFIM 
systems coupled by a transverse coupling t±_ (or Ji_). The T = 
Ordered Phase of the uncoupled TFIM is now extended, with 
a phase boundary which has a value Tc for the model at <; = 
and a T = quantum critical point (QCP) aX g = Qc- The 
hashed region immediately to the right of the ordered phase 
and just above the QCP is the gapless quantum critical re- 
gion (described in the text). The dashed line represents a finite 
T crossover from the quantum critical region to a Disordered 
phase. 

In keeping with the fact that the passage to higher 
dimensions is essentially a nonperturbative phenomenon 
[7], we use the RPA method to study the physics of di- 
mensional crossover. Put another way, we are interested 
in gaining an understanding of how lower dimensional sys- 
tems (here, our one-dimensional strongly correlated chains), 
when coupled to one another, go from being nearly iso- 
lated to a anisotropic strongly coupled system in higher 
dimensions. It is worth noting that, while the mean-field- 
like approach of RPA is exact only in infinite dimensions 
(i.e., infinite coordination number), its application to the 
physics of coupled quasi-lD spin systems for small coordi- 
nation numbers (i.e., lower dimensions) has met with con- 
siderable success Unmi]. Further, while a naive mean-field 
treatment of single particle-hopping between fermionic chains 
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is not possible as a single fermion operator has no well- 
defined classical limit [Tj , we are able to treat two-particle 
hopping processes between our underlying chain and lad- 
der systems via RPA by working with an effective theory 
in terms of pseudospins. This is justified because the pres- 
ence of the large on-site Hubbard coupling in our under- 
lying model makes the single particle hopping irrelevant 
(in a RG sense), while two-particle processes (including 
hopping terms) can be crucial in determining the phase 
diagram Pfl . A full treatment including single-particle hop- 
ping will require a treatment using chain-dynamical mean 
field theory (c-DMFT) [42] and wiU be the focus of a fu- 
ture work. 

For a transverse coupling 0±, the RPA method in- 
volves computing the dynamical spin susceptibility x of 
the coupled system in the disordered phase as 



x{uj,k,k^) = [x^^{ij,k)-o^{k^)Y 



(29) 



in terms of the frequency uj, the longitudinal and trans- 
verse wave-vectors k and fej^ respectively. xiD is the dy- 
namical spin susceptibility of a single TFIM system, to be 
calculated assuming incipient order along the direction 
in pseudospin space 



poo 

/ dte<-'-'^^\[T^{t,n),T%0,0)]) 



(30) 

and the transverse coupling 0±{k±) ^ z±0±{k± = 0), 
for each TFIM system having a coordination number of 
z± . Then, a divergence in the dynamical pseudospin corre- 
lation function x{^^ k, k±) signifies an instability towards 
the formation of an ordered state. In this way, it is pos- 
sible to compute the quantities aX T = 0, Tc for the 
case of 5 = Z\/(y — J/4 — 2P) — (single pseudospin 
chain mass A = \V - J / A - 2P - 2t\ for the TFIM), the 
shape of the phase boundary near the QCP, the dynam- 
ical spin susceptibility x(w, fc, k±) at the QCP as well as 
the dispersion in the transverse directions for small 
and close to the QCP pl5]. We focus our attention mainly 
on, and in the neighbourhood of, the QCP in order to 
assess the role played by the critical quantum fluctua- 
tions in determining the physics of deconfinement and di- 
mensional crossover in our system. In this, we will often 
be aided by the integrability and conformal invariance of 
the unedrlying pseudospin model (e.g., for the TFIM, at 
2t = V-J/A- 2P, O±=0 (SUEZ]). We can then exploit 
the nonperturbative results for a single pseudospin system 
(as long as the mass scale A << 1), while dealing with the 
physics of the transverse couplings at a mean-field level. 

We learn from the results presented below that, in 
passing from the interior of the ordered phase towards 
the phase boundary, the excitation gaps decrease together 
with a gradual growth of the dispersion in the transverse 
directions. The spectrum is gapless at the QCP. At low 
T, directly in the region above the QCP, the spectrum 
and dynamics are mainly governed by the QCP while the 
thermal excitations are described by the associated contin- 
uum quantum field theory [57] . The dimensional crossover 
is thus characterised by the growth of the dispersion in the 



transverse directions close to the QCP, while the T = 
deconfinement transition is characterised by the vanishing 
of all mass gaps at the QCP. 

We can, thus, begin with the case of transverse cou- 
plings t± and J± in a system of many TFIM chains (as 
given in cqun. lj^T]) ). As a detailed study of a generic sit- 
uation of coupled TFIM systems has been carried out in 
Rcf . f [191135] ) . it is sufficient to quote the main results here 
for the case of t±. As discussed below, this is done since 
the results for the J± coupling are found to be qualita- 
tively similar. The critical temperature Tc for the case of 
g ~ is obtained from the Xid{^ = 0, A; = 0) at finite T 



Xid(c^ =0= k)= 



C2 



2t:T 



N-7/4 



{V - J/4 - 2P) ' {V - J/4 - 2P) ' 



(31) 

where the constant C2 = sin(7r/8)B^(l/16, 7/8) and B(a;, y) 
is the Euler Beta function. Thus, we find, from the condi- 
tion for the divergence of the susceptibility x 



iy - J 1 4 - 2P) 



C3( 



iy - J/ A - 2P) 



4/7 



(32) 



where the constant C3 = C2/{2tt) — 2.12. 

We next compute the critical value of the transverse 
coupling gc aX T = 0. For this, we can use the expression 
for slightly off-critical xiD for the case when the mass of 
the single TFIM system is very small (m — Ac « 1). 
This, for small w, is given by 



XiD{i^,k) 



^2 _ („^)2 _ ^2 



(33) 



where the velocity v — {V — J / A—2P)a (where a is the lat- 
tice spacing) and Zq = 1.8437. Then, from the divergence 
of the susceptibility of the coupled system, x('-^i k^ kj_) 



Xioi^^k) = z±_t^{k±_ = 0), 



we get, for the case oi lo = {) = k, 



ly - j/4 - 2P) 



,7/4 



Z\t_ 



(34) 



(35) 



where we've dropped the argument of k^ = in t_L for 
the sake of convenience. Thus, we get 



5c « ci ( 



z±_t^ 



{V - J /A - 2P) 



)4/7 



(36) 



4/7 

where the constant ci ^ Zq = 1.42. Precisely the same 
expression for the mass in the ordered phase and very 
close to the QCP, A = gc{V - J/A - 2P), is also obtained 
by carrying out a self-consistent treatment of the effective 
magnetic field, h = z±t±{T^) in the TFIM for a single 
chain. For this, one uses the slightly off-critical suscepti- 
bihty for the ID TFIM given earlier (eqn.Q) with the 
mass A replaced by A{1 + {h/{V - J/A - 2P)f) [43]. 
From this result, the authors of Ref.([Tn]) concluded that 
the dispersion in the transverse directions in the ordered 
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phase and close to the QCP is much stronger than that 
deep in the ordered phase. 

The susceptibility x for the coupled system at the QCP 
for small fcj^ can also be computed by using the relation 
for the slightly off-critical xiD given earlier, eqn. (P5|) . to- 
gether with the relation Al = gl^'^iV - J/4 - 2P)t±_ in 
eqn.jSll). This leads to 



1/4 



Zo{V-J/A-2P)g, 
uj^ - {vk f - {vi_ ■ fe_L)2' 



(37) 



where \v^_\^ = {Zo{V-J/4.~2P)gl''^ /2)dH^{k^ = {))/dk\ 
is gained by a Taylor expansion to second order 19,35 . 
The shape of the phase boundary at low T can now be 
determined by using the xiD of the TFIM at low T [27] 



(38) 



where — ^e^/"^ is the dephasing time due to quantum 
fluctuations. Then, for = = fc in Xi-Di the ean. (p4)l 
gives 



In T = In m 

T 



Inyl, 



(39) 



Z(\t I 



where A = f (^^y/iy 
above has an approximate solution [191135 



1)^/^. The expression ([39]) given 



(40) 



ln(l/yl)-lnln(l/yl) 



This relation gives us the shape of the phase boundary for 
low T and close to the QCP. We have, therefore, derived 
the important features of the T — g phase diagram for the 
case of the transverse coupling (Fig.([l|)) given above. A 
RPA calculation for the other transverse coupling, J^, for 
critical TFIM chains (i.e., the Hamiltonian ([^Tjl with V ~ 
J/4 — 2P = 2i, U±^ — J± and t± — 0) can also be carried 
out ^9\. This is guaranteed by the integrability of this 
Hamiltonian [44] and the fact that it falls into the same 
universality as the TFIM. The RPA calculation thus leads 
to an similar set of relations for gc, Tc, the susceptibility 
X, dispersion in the transverse directions as well as the 
shape of the phase boundary close to the QCP to those 
obtained earlier, but with t± replaced by J± everywhere. 
In this way, we obtain essentially the same T — g phase 
diagram for the case of the J^ transverse coupling; the 
only notable difference is that the spectrum of the ordered 
phase is now obtained from the exact solution of the ID 
TFIM in a longitudinal field ,,44] . 

We now turn to a discussion of the results of a RPA 
treatment for the case of coupled 2-leg TFIM ladders, 
by using the LEEHs obtained in various limits in sec- 
tion III. For the case of the strongly coupled 2-leg Lad- 
der LEEH with the coupling Jj^ > (equ.((25])), as we 
again find the TFIM Hamiltonian as the effective the- 
ory, we can safely conclude that an RPA treatment for 
many such coupled ladders will lead to the same results 
as those given above and a phase diagram identical to 
Fig.®. A RPA treatment of the LEEH derived for the 



case of J^ < fequn. ([^5)) ) can be carried out easily for 
some special cases: ii)h = = t_L, (ii) t± = {J^ / J±){J'j^ — 
2/i2)/(j2 - 4/i2) and (iii) J^ = 2h^. We discuss each of 
these in turn. 

For t± = = h, the LEEH for J± < reduces simply 
to that of the ID Ising model with an exchange coupling 
—J^/2J±. We can now treat interchain coupling terms 
^iij.Jjj, + ^la^lb) ^^'^y^la^lb (where a, b are chain in- 
dices, i ^ J^IJ\ and U ^ J^/Jl) in turn via RPA. In 
this program, we can assume order along a given direction 
in pseudospin space and then replace the appropriate in- 
terchain coupling term by an effective field term (in the 
spirit of a mean-field treatment) and solve the Hamilto- 
nian in a self-consistent manner. Thus, we can see that for 
the t spin flip term, such a mean-field treatment assumes 
order along say, and thus has an effective magnetic 
field h = z±t{£^^) (where z± is the coordination number 
of any ladder system) . This leads to the self-consistent ef- 
fective Hamiltonian of the coupled ladder problem again 
taking the form of the ID TFIM. Computing the critical 
value of the gap at the QCP in a self-consistent manner 
(as discussed earlier) gives gc = /{2J±){2z±tJ±/ J^)^^'^ . 
In this way, it is clear that the self-consistent solution of 
this effective theory will again give rise to a phase diagram 
like Fig.®. 

An identical calculation for the U Ising transverse cou- 
pling term can be carried out for the case of i_L = ( J^/ J_l)( J^ — 
2/i2)/(j2 - 4^.2). This is the case of the critical TFIM in 
a longitudinal field, which is integrable and belongs to 
the same universality class as the TFIM [44 . The RPA 
calculation assumes order along and has an effective 
(self-consistent) magnetic field h = z±U{£,^) in a ID Ising 
model in a longitudinal field. Using the exact solution of 
the problem [44] , together with the divergence condition in 
the RPA, gives us the critical gap as gc = J^/(2 J_l)(2z^C7Jj^/ J^)*/'^. 
Again, we reach qualitatively similar conclusions with re- 
gards to the phase diagram for this transverse coupling. 
Finally, we note that the LEEHs for both the cases of 
J_L < 0, J2 = 2/i2 and tj_ » J, J± are those of the XY 
model in a transverse magnetic field; a detailed RPA in- 
vestigation of this problem will be presented elsewhere. 



5 Comparison with recent numerical works 

We present here a discussion of the relevance of our work 
to some numerical investigations that have been carried 
out on strongly correlated single chain, 2-leg ladder sys- 
tems as well as coupled TFIM systems. We begin with 
a discussion of the early work of Capponi et al. 8J. The 
study assessed the effects of long-range Coulomb interac- 
tions on the phase diagram of a finite-size one-dimensional 
system of spinless electrons at 1 /2-filling. The authors con- 
cluded from an exact diagonalisation analysis that for in- 
termediate strengths, the presence of extended range in- 
teractions caused an enhancement of the metallic nature 
of the system; this is in contrast with the fact that in 
the thermodynamic limit, such a system would always be 
driven by the logarithmic divergence of the long wave- 
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length part of the interactions towards insulating charac- 
ter. While the metallic phase was observed to have a van- 
ishing gap, it did not agree with the predictions of confor- 
mal field theory. Increasing the strength of the Coulomb 
interactions caused a crossover towards a localised 2kF 
CDW phase. Note that a gapless metallic phase arising 
from a Wigner CO phase as the strength of the Coulomb 
interactions is reduced is not in contradiction with our re- 
sults for the charge sector of the single chain in the strong- 
coupling regime: the effective TFIM model at T = has 
a gapless QCP emerging from a CO phase as the nearest 
neighbour hopping t grows to a critical value. Capponi et 
al., however, found no signatures of the Peierls CO phase 
for the single chain discussed in section II. 

There are a few notable works on quasi-lD strongly 
correlated models at 1 /4- filling with extended interactions. 
Riera et al. [T^ studied the case of 1/4- filled chain/ladder 
Hubbard and t — J systems, but which also include Hol- 
stein and/or Peierls- type couplings to the underlying lat- 
tice. Their findings reveal coexisting charge and spin or- 
ders in both chain and ladder systems. Specifically, for the 
case of their chain system, by keeping only on-site and nn 
repulsion together with an on-site Holstein-type coupling 
of the electronic density to a phonon field, their phase di- 
agram (Figs. 4) reveal separate phases with Wigner as well 
Peierls- type charge order. This is in keeping with our find- 
ings, but the origin of the Peierls order in the two cases are 
different: in our work, it originates from the competition 
of the nnn coupling V2 with the nn coupling Vi , while in 
their work, it needs the Holstein coupling to the lattice. 
Riera et al. find a similar phase diagram (Figs. 5) for the 
case of an extended t— J model (i.e., including nnn t and 
J couplings) with a Holstein coupling. The addition of 
a Peierls-type coupling leads to a spin-Peierls instability, 
i.e., the formation of a dimerised spin order, which coex- 
ists with the Peierls-type charge order. Again, while this 
matches our findings, the origins are different. Qualita- 
tively similar conclusions are also reached by the authors 
in their study of an anisotropic 2-leg t — J ladder with an 
extended on-chain nn coupling and Holstein/Peierls-type 
lattice couplings (Figs. 2 and 3). 

Vojta et al. [45 studied the problem of a strongly 
correlated electronic problem at 1/4- filling and with ex- 
tended Hubbard interactions (keeping only a nn repul- 
sion) using the DMRG method. The phase diagram they 
obtained contains several phases with charge and/or spin 
excitation gaps. While a comparison of our work with this 
study is hindered by the fact that the DMRG analysis 
does not have the crucial element of the nnn repulsion 
(V2 in our work). Fig. 2 of that work reveals that for the 
case of [/ » Vi > t, the authors indeed find a charge- 
ordered CDW state (i.e, the Wigner charge ordered state 
of discussed in section II) with an excitation gap in the 
spin sector as well. This is in conformity with our find- 
ing of a Wigner charge-ordered state with a spin gap for 
the case of the Up coupling being the most relevant un- 
der RG. Further, the t± of that study corresponds to the 
single-particle hopping between the legs while our work 
has focused on the effects of two-particle hopping. Finally, 



with the on-site Hubbard coupling, U, being the largest 
in the problem, we are unable to see any phase-separated 
state in our phase diagram (as observed by Vojta et al.). 

We end by commenting on a very recent DMRG stud- 
ies of Konik et al. [1^ on coupled TFIM systems in an 
effort at studying two dimensional coupled arrays of one 
dimensional systems. By starting with a reliable spectrum 
truncation procedure for a single TFIM chain (which re- 
lies on the underlying continuum ID theory being either 
conformally invariant or gapped but integrable), the au- 
thors then implement an improvement of their DMRG 
algorithm using first-order perturbative RG arguments. 
Their results for a J± coupling of the TFIM chains con- 
firms the accuracy of the RPA analysis of Ref.f[19|) and 
the present work in computing quantities like the single 
chain excitation gap (which is found to vanish at a criti- 
cal J±) and dispersion of excitations in the coupled sys- 
tem as a function of J±. Their results indicate that the 
RPA method and the DMRG analysis agree very well upto 
values of J± of the order of the gap. This method also ap- 
pears to give accurate values of critical exponents related 
to the ordering transition. Thus, this numerical approach 
appears to provide a confirmation of the interplay of the 
QCP in the TFIM and the transverse coupling in driving 
the dimensional crossover and deconfinement transition. 
Such an approach, therefore, holds much promise for the 
numerical investigations of such phenomena in systems 
with similar ingredients. 

6 Conclusions 

To conclude, we have explored the strong-coupling limit 
of strongly correlated 1/4- filled single chain and two- leg 
ladder models using a variety of methods. The charge sec- 
tor of a 1/4-fillcd ID model of electrons with extended 
interactions is, in the regime of the on-site Hubbard term 
being the largest and the hopping strength the smallest, 
found to lead to an effective theory of the ID transverse 
field Ising model. The ordered phases in the charge sector 
are found to belong to either the Wigner or Peierls type 
CO. The two kinds of CO are then found to give rise to 
a gapless AF ground state and a dimerised state respec- 
tively in the spin sector. The integrability of the ID TFIM 
allows us to make considerable progress in computing var- 
ious thermodynamic quantities (e.g., response functions), 
especially in the quantum critical region lying just above 
the T = quantum critical point (QCP). 

Strongly coupled 2-leg ladder systems composed of such 
chains are also studied in various limits, with the low 
energy effective theory generically being found to be de- 
scribed once again by an exactly solvable ID model with 
a QCP. The varying of interchain couplings in the 2-leg 
ladder is found to give rise to a variety of charge ordered 
phases (e.g., in chain Wigner CO, rung-dimer as well as 
orbital antiferromagnet type charge order) . This is in con- 
formity with our recent studies on the weak coupling phase 
diagram for a coupled 2-leg TFIM system ^35j. Doping 
such ladders with holes is also found to give rise to su- 
perconductivity. RPA studies on the effects of transverse 
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couplings connecting many such chain and ladder systems 
presented a generic phase diagram for the coupled sys- 
tem containing an ordered phase extending to finite T 
and with a phase boundary ending in a quantum critical 
point. These calculations also stressed the importance of 
the role of the QCP in the mechanism responsible for the 
dimensional crossover in the quantum critical region ly- 
ing at finite T just above the QCP and the accompanying 
T = deconfincnicnt transition. 

Significantly, the favourable comparison of our find- 
ings, as discussed earlier, for the strongly coupled sin- 
gle chain and 2-lcg ladder systems with experimentally 
observed phenomena in prototype examples like organ- 
ics (TMTSF) and SruCuuOn respectively strongly sug- 
gests that these systems lie in close proximity to an under- 
lying QCP associated with charge order. This constitutes 
a significant advance in our understanding of the physi- 
cal responses of these systems in a new theoretical frame- 
work. Further, the robustness of the dimensional crossover 
mechanism found in this work leads us to conclude that 
critical quantum fluctuations associated with a QCP can 
quite generically enhance the dispersion in the transverse 
dimensions for anisotropic systems, i.e., facilitate the pas- 
sage from the gapped phases of lower dimensional systems 
to the gapless phases of the coupled system in higher di- 
mensions. 

We thank E. Miiller-Hartmann, G. I. Japaridze, F. Pranchini, 
L. Dall'Asta, S. Basu, S. T. Carr, A. Nersesyan and M. Fabrizio 
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